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A FORMULA FOR THE RELIABILITY OF A d-DIMENSIONAL 
CONSECUTIVE-fc-OUT-OF-n:F SYSTEM 

SIMON COWELL 


Abstract. We derive a formula for the reliability of a d-dimensional consecutive- 
fe-out-of-n:F system. That is, a formula for the probability that an ni X ... X 
array whose entries are (independently of each other) 0 with probability p and 
1 with probability g = 1 —p does not include a contiguous si X ... X subarray 
whose every entry is 1. 


1. Introduction 

Consider a pipeline punctuated by several pumping stations. Suppose the pump¬ 
ing stations are redundant in that, if one of them fails, then its predecessor will 
be strong enough to pump the fluid past it to the next pumping station. Perhaps 
each pump will even be able to compensate for the failure of its successor and its 
successor’s successor, but perhaps not for the failure of the first 3 pumping stations 
in the chain of its succesors. In this case, the system fails if and only if among 
the sequence of pumping stations there exists a consecutive run of 3 or more failed 
pumping stations. This is an example of what is known as a consecutive-A:-out-n:F 
system. 

Suppose that such a system having n nodes fails if and only if k consecutive 
nodes fail, and suppose that any given node in the system works correctly with 
probability p € [0,1], independently of the other nodes. Then the probability that 
any given node fails is g = 1 — p. The Reliability of the system is the probability 
R{k,n\q) that the system does not fail. Then R{k,n',q) = 1 — P{k,n',q), where 
P{k,n]q) is the probability that the sequence of n nodes includes a contiguous 
interval of k or more failed nodes. Here we are using notation as in [DB2I4) . 

As noted in [DB214j . the concept of the reliability of a consecutive system was 
introduced to Engineering by Kontoleon in 1980 |K80) . In the following year Chiang 
and Niu [CN8I| discussed some applications of consecutive systems. The concept 
has been generalised in several directions, for instance to systems deemed to have 
failed if and only if they include: k consecutive failed components or f failed 
components [T82| ; k consecutive components of which at least r have failed |G86| ; 
at least m non-overlapping runs of k consecutive failed components |P90) . In 2001 
Chang, Cui and Hwang published a book on the subject [CCHOl) . Four reviews 
are also available [CFK95] . |CH03] . |E10] and |T15) . 

An exact formula for R{k, n; q) first appeared in [D1738| . De Moivre’s approach 
depends on deriving the generating function (see [W2005] l for P{k, n; q). This work 
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has been rephrased in a modern style in [U37] . Much later a different exact formula 
for R{k,n;q) using Markov chains was found |F86) . [CF89) and [FH87] . 

Both of these exact formulae for R{k,n\q) are difficult to use in Engineering, 
where values of k and n may be very large. Efficient algorithms have been found 
1 [CKS08| and |HW95j l but even with modern computing resources, evaluating the 
formulae takes a lot of time and memory. This problem has been addressed at length 
in the study of upper and lower bounds for n; q). Such bounds have been found 
which are close approximations to the exact value of i?(fc, n; q) and which are also 
much easier to evaluate than the exact value oi R{k,n\q). In |DB114j . |BD14j and 
[DB214j the authors compare many such results from the 1980s to the present day. 

Dau§ and Beiu were the first to use De Moivre’s original formula to derive upper 
and lower bounds for R{k,n;q). As shown in [DB214) this method compares very 
favourably with the previous upper and lower bounds, especially for small g, and 
also with the exact algorithms given in |CKS08) and |HW95) (using their bounds 
is about 100 times faster than the algorithm from [CKS08| . and 1000 times faster 
than the algorithm from |HW95] b 

One possible extension to higher dimensions is as follows: Given d S N and 
ni, n 2 ,..., G N, suppose that some system is represented by an ni x n 2 x • • • x nd 
array of Os and Is. Suppose that each node in the array is 0 with probability 
p G [0,1] and is 1 with probability q = 1 — p, independently of the other nodes. 

Suppose that si, S 2 , ■ ■ ■, Sd G N are given, with Sr < Ur for each r, and suppose 
that the system is deemed to have failed if and only if the d-dimensional array 
includes a contiguous si x S 2 x • • • x Sd subarray of Is. Define the Reliability of 
the system as the probability R{si,..., Sd,ni,..., Ud] q) that the system does not 
fail, and let P(si,..., Sd, «i, ■ ■ ■, g) = 1 - P(si,..., Sd, ui,..., Ud; g), so that 
P(si,..., Sd, 111 , ■ • •, rid', q) is the probability that the array includes a contiguous 
Si X S 2 X ■ • • X Sd subarray of Is. 

Let us survey in chronological order some recent papers dealing with the higher 
dimensional case. In [SL90j the 2-dimensional special case i?(s, s, n, n; g) is treated. 

Upper and lower bounds are found by reducing to the 1-dimensional case. Using 
empirical approximations to the exact value, i?(s, s, n, n; g) is plotted as a function 
of g for n = 5 and s = 2,3 and 4, for n = 10 and s = 2,4,6 and 8, and for 
n = 50 and s = 2, 5,10 and 20. The authors do not attempt to find an explicit 
expression for R{s,s,n,n',q). In |BKP92j the 2-dimensional case is considered, for 
R{si, S 2 ,ni,n 2 ', q) with (si,S 2 ) G {l,ni} x { 1 , 712 }. Exact formulae are derived by 
comparison with the 1-dimensional case. Similar results are also given for a cylin¬ 
drical version of the 2-dimensional array. The 2-dimensional case R{s, s,n,n', q) of 
[SL90| is revisited in |KPP93j . Again, upper and lower bounds are given, but in the 
more general case where the probabilities of distinct nodes of the system failing are 
not necessarily equal. For the general 2-dimensional case i?(si, S 2 , ni, 712 ; g), recur¬ 
sive formulae for the exact reliability have been given by Yamamoto and Miyakawa 
[YM95| and also by Noguchi, Sasaki, Yanagi and Yuge [NSYY97] . The result pre¬ 
sented herein is distinct from their work, in that we give a closed-form formula for 
the exact reliability as a polynomial in q, for a general system, of arbitrary dimen¬ 
sion. The 3-dimensional case is treated in [GESNIO] . where exact formulae for the 
reliability are given in special cases analogous to those studied in |BKP92j . 

In the present article we derive general closed-form exact formulae for R{si,..., Sd,ni,..., Ud', q) 
and for P(si,..., Sd, ni,..., Ud', g), where the dimension d G N is arbitrary. As far 
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as the author knows, no such formulae were known before for d > 1, even in the 
case d = 2. Indeed, in [KPP93] the authors write “It is very difficult (probably 
impossible) to derive simple explicit formulas for the reliability of a general 2D- 
consecutive-/c-out-n:F system.” 


2. Results 

Theorem 2.1 (Main result). Let the terms 

d, 

^li ' ■ • ^d-i 
Til j ■ • ■ 


R{si,...,Sd,ni,...,nd-,q) 
P{si,...,Sd,ni,...,nd;q) 
be as in the introduction, and let 

d 

FI = J^{1,... ,nr - Sr + !}• 

r—1 


Then 


(1) R{si,...,Sd,ni,...,nd-,q) = 


1+ E (-1)'"' n 

JGViE)\ld J'&ViJ)\ld 


(—1)1'^ I In max ^0, Sr 


and 


( maxCr 
eeJ' 


min 

e^J' 



(2) P(si,...,Sd,ni,...,nd;q) = 


- E (-1)'"' n ®^p 

Je'P(E)\0 J'G-P(J)\0 


(-l)l'^'lln 


max I 0, Sr — I 


maxcr 

\eeJ' 


— min 
eG J' 



where 'P{A) denotes the power set of the set A, that is, the set of subsets of A, and 
0 denotes the empty set. 


Theorem 12.II implies the following combinatorial result: 


Corollary 2.2. With all terms as in Theorem \2.1\. the number a{si ,..., s^, ni,..., Ud) 
of failed systems among the 2 Ar=i'^r ppggiyig gygtems is 


q(si, • ■ •, Sd, n\, • ■ •, ndf) 
\\J\ 


_ 2 nr=i”'- ^ (—1)''^' exp 

J&ViE)\li J'G-P(J)\0 


(_l)l'^ I C2) TT max ( 0, Sr — ( maxCr — mine, 
V V eGJ' eGJ' 

r=l ^ ^ 


Proof. In case p = q = ^, all possible systems occur with equal probability, so the 
probability measure on the power set of the set of all systems becomes the counting 
measure, divided by a factor of 2^^=^ □ 
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Setting d = \, ni = n and si = 2, 3 and 4 in Corollary 12.21 yields the sequences 
[SKW99) . [W199] and [W299) . respectiyely. Setting d = 2, ni = 712 = n and 
Si = S 2 = 2 in Corollary 12.21 yields the sequence |C15] . whose complementary 
sequence is [H08]. 

Proof of Theorem \2.1i Call a system e an “elementary failure case” if and only if 
it includes a contiguous si x ... x Sd subarray of Is, and its other elements are all 
0. We identify the set of all elementary failure cases e with the set E, according to 
the rule 

e yA (ei, ...,ed) 

if and only if the (ii,..., id)-th element of e is 1 precisely when, for all r € {1 ,..., d}, 

67 ^ ^ hj- ^ Cy “t“ Sy 1 . 

For example, in the two-dimensional case, we identify an elementary failure case 
with the indices (ei, 62 ) of the top left corner of its contiguous subrectangle of Is. 
Denote by G the simple acyclic directed graph whose vertices are all the possible 
systems, and in which system & is a direct successor of system a if and only if b is 
obtained from a by changing a single 0 element to 1. Then the failed systems / 
in G are precisely those vertices having some elementary failure case e G £1 as an 
ancestor (including the possibility that / = e.) Therefore if F denotes the set of 
all failed systems / (in terms of probability, F is the event that the system fails,) 
then 

F = (^ the set of descendents of e, 

eeE 

where by “b is a descendant of a” we mean to include the possibility that b = a. 
Then 

P(si, ...,Sd,ni,...,nd;q)= P{F) 

= P I LJ the set of descendents of e 
VeeE 

= E (- 1 )'-"^'^ n the set of descendents of e 
JeV{E)\9 VeeJ 

by the Inclusion-Exclusion Principle for the probability measure on the power set 
of the set of all possible systems. But 

Pi the set of descendents of e = the set of descendents of hcd(J), 

egj 

where by hcd(iF) we mean the highest common descendant of any given set K 
of vertices in G. That is, in the partial order on the vertices of G defined by 
a < 6 precisely when a is an ancestor of &, hcd(£f) is the least upper bound of K. 
Moreover, 

hcd(J) = [J e, 

eeJ 

where by the union of several ni x • ■ • x systems we mean the system whose 
generic element is 1 if and only if the corresponding element in at least one of those 
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systems is also 1. Thus 

( 3 ) 


Jev{E)\ii 

Jev{E)\ii 



However, for any vertex a of G, 

P(the set of descendents of a) 


E 

6 is a descendent of a 


E 


nin2...n^-k k 
y *1 5 


6 is a descendent of a 


where k = k(b) is the number of failed nodes in the system b. Gathering like powers 
of q in the above sum, we have 


nin2...nd / 

E 


P(the set of descendents of a) = 


k—k{a) 

nin2...nd 

^k{a) 


nin2---nd 

k — k{a) 


= r- E 

k—k{c 

711712.■■rid —k{a) 

E 

3=0 

— qk{a) _|_ ^■^nin2...nd-k{a 

= q^^°‘\ 


711712 .. .Ud — k{a) 
k — k{a) 


pnin2...nd-k ^k-k{a) 


= q^^ 


ni7l2 ■ ■ .nd ^(®)^^nin2...nd-fe(a)- 


That is, for any vertex a of G, 

P(the set of descendents of a) = q^‘^°‘\ 

where k{a) denotes the number of failed nodes in the system a. Therefore from ([3]) 
we obtain 


( 4 ) 

P(si,...sd,ni,...7id;g) 


E 

Je'P(E)\0 


(_l)PI+ip gg^ Qf descendents of 


= ^ (_l)PI + lgfe(Ueeje)^ 

JeP(E)\0 


The next task is to compute the number k (IJeeJ ®) system IJggj e, 

where J is some nonempty subset of the set E of elementary failure cases e. We 
use the Inclusion-Exclusion Principle again, this time for the counting measure on 
the power set of the set of elements of a system, to obtain 

(5) ^ f u = E f n ■ 

VeGJ / J'eT’(J)\0 VeeJ' / 
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Each elementary failure case consists entirely of Os except for a contiguous si x 
... X Sd subarray of Is. Therefore PleG J' ® consists entirely of Os except for a 
(possibly empty) contiguous subarray of dimensions ti x • • • x td, where for each 
r e 


tr = max 



I min (cr + s 
\eeJ' 


r 


1 ) 


— max Cr + 1 

e^J' 


= max 



( max Cr 
eGJ' 


— min e 

e£j' 



The volume of that contiguous ti x ■ ■ ■ x td subarray is equal to k (HeGJ' ®)’ 
is, 



Considering (|4]), ([5]) and (|6]), we have 


( maxer 
eeJ' 


— min 

e^J' 



P{si,...,Sd,ni,...,nd;q) = ^ (_l)l‘^l+ig'=(Ueer«) 

JeP(E)\0 


E (-l)'^'+'exp 

Hq) E 

[n») 

JeP{E)\9i 

J'e'P(J)\0 

VeeJ' / 


= 51 (-l)'^'+'exp 

JeViE)\<!i 


ln(g) 


E (-1 

J'GT’(J)\0 


)|T|+i J^niax To, Sr - ( 


max tr 

\e^J' 


— min 
ee J' 



= - E (-1)'"' n 

J&V{E)\<Ii J'eV{J)\<!i 


(-i)Vii„(i)n, 

\ *1 / j^—i 


max I 0, s 


- 


max Cr — mm 

\eGJ' eGJ' 


which proves equation ([2]). Equation m follows immediately. 


□ 


As an example, we compute the reliabilities of 3 consecutive systems, having 
dimensions 1, 2 and 3: 


^(l,2;g) = 1 - 2g + 9^ 
i?(l, 2, 2, 3; <7) = 1 - 4q2 + 2q^ + V - + 9®, 

i?(l, 2, 3, 2, 3,4; q) = 1 - + 4g® + 4q^ + 4q^° - 8g“ + ISq^'^ - 16q'^'^ - 16q^^ 

- 12gi® + 409^^ + 4gi® - 8q^^ - 8q^° - 12q^^ + 20q'^^ - 8q^^ + q^^. 


See also Figure 1. 
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Figure 1 



R(\,2-,q) 

- R(l,2, 2, 3;g) 

i?(l,2, 3, 2, 3,4; q) 


3. Conclusion 

We have provided a novel formula for the reliability of a general d-dimensional 
consecutive-A:-out-of-n:F system, as an exact polynomial in q. We believe ours 
to be the first general exact closed-form formula published for the reliability in d 
dimensions. This answers an open question in Engineering. 
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